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ABSTRACT 

This thesis demonstrates the significant benefits in the 
lynamic stability improvement of a power system, consisting of a 
iynchronous generator connected to an infinite bus through a 
Lransmission line, by reactive power modulation at the generator 
terminal and at the mid— point of the transmission line, 
iimul taneously . Reactive power modulation at the generator 
terminal is achieved by incorporating power system stabilizers 
tfhile at the mid point of the transmission line it is achieved 
>y means of Static Var compensator having an auxiliary 
:ontroller. This thesis proposes a method to design and 
inordinate P.S.S. and S.V.C. stabilizer. The benefit derived in 
Lhe robustness by such coordinated application of PSS and SVC is 
also demonstrated. 

A suitable model for dynamic stability analysis of a single 
nachine connected to the infinite bus through a transmission line 
tfith Static Var compensator at the mid point of the transmission 
line is also proposed. 

The procedure followed here is simple. A set of P.S.S. is 
designed assuming that S.V.C. is not having an auxiliary 
controller. Similarly a set of S.V.C. stabilizer is designed 



xiv 


assuming that PSS is not connected in the system. The design 
point chosen is the same for both designs. The PSS and SVC 
stabilizer so designed are then coordinated so that when both are 
connected in the system a further improvement in the region of 
robustness is obtained. Pole placement technique is used for the 
stabilizer design and a sector criterion for eigenvalues is used 
for roubustness. 

The proposed method has been tested for various 
combinations of control signals- 



CHAPTER 1 


INTRODUCTION 

1.1 G«nAral 

Power system engineers have to operate their systems in such 
a way that they can meet the load demand on it in the most 
economical manner, ensure good quality and high reliability 
supply and meet various contigencies in operation effectively. 
This is a difficult task. It has been compounded by the fact 
that power systems have grown rapidly and have become very 
complex, to meet with the ever increasing load demand. 

Enhancement of power system stability is a major factor in 
improving the reliability and quality of the system. Power 
system engineers speak of three types of stability. The firstis 
steady state stability which puts a limit on the maximum power 
transfer capability of the system. The second is transient 
stability. It is the capacity of the system to withstand major 
shocks like loss of a generating unit, loss of a major 
transmission line or a severe network fault without losing 
synchronism. And the third is dynamic stability. We confine 
ourselves in this thesis to the problem of dynamic stability. 

In a power system small magnitude random changes in load 
take place at all times, with the accompanying adjustment of 
generation. When load on the system changes, the generators in 
the system move from one operating point to another. The 
transition following a system perturbation is oscillatory in 
nature, but if the system is stable, the oscillations will be 
damped out and the system settles down at a new steady state. 



The ability of all the machines in the system to adjust to small 
load changes is called dynamic stability of the system. If 
damping is insufficient, the oscillations may grow and the 
growing oscillations on a heavily loaded transmission line impose 
limitations on power transfer limit, leading to the tripping of 
the transmission line or forcing the under utilization of the 
transmission line capacity. Low frequency oscillations of 
sustained or growing nature have been reported in the literature 
by the power utilities E1,2D. These oscillations typically 
occur in the low frequency range of about 0.2 to 2.0 Hz. These 
oscillations in power can be related to oscillations of rotors of 
synchronous machines in the system relative to one another. In 
the literature these low frequency oscillations are called by 
different names like zero mode, ih,^ltjbial mode, elect r omechanical 
mode or rotor mode. 


1.2 EffActs of Thyrlst/Or typA Exclbatlon Sys-Lam on Dynamio 
Stability 

The modern practice is to use fast acting, high gain, static 
excitation systems for synchronous generators. They improve the 
transient stability limit of the power system, but have 
•deleterios effect on dynamicstabil ity, resulting, often, in 
lightly damped and even sustained oscillations. This is more 


pronounced during heavy loading conditions. 

de Mellow and concordia in their celebrated paper C3D have 
provided an insight into the effects of thyristor type excitation 


systems and established understanding of 


stabi 1 izing 


requirements of such systems. 


the 



1.3 Methods of Damping Low Frequency Oscillations 

Power system engineers started serious research towards 
finding out means to damp out low frequency oscillations after 
the reseach work reported by de Mellow and Concordia C33. As 
already mentioned, thyristor type excitation system introduces 
negative damping in the system. To offset this effect and to 
improve system damping in general, artificial means of producing 
torque in phase with speed are introduced using supplementary 
signals in the excitation control system of generators. 

Additional networks are used in the excitation system to 
generate stabilizing signals. They are called ’Power System 
Stabilizers’ (P.S.S.). P.S.S. introduce stabilizing signals in 
the excitation system. There is always a strong coupling between 
the voltage magnitude and reactive power in a power system. 
Hence P.S.S. essentially modulate reactive power to improve 
dynamic stability of the system. 

Another method of damping the low frequency power 
oscillations by modulating reactive power is gaining momentum in 
recent years. Recent advances in thyristor technology have 
resulted in the emergence of a new class of devices known as 
’Static VAR Compensators' (S.V.C.). The main application of such 
devices is to improve system voltage profile on the system by 
providing reactive power compensation. 

Though the main function of S.V.C. is to control voltage, it 
has been found that a significant contribution to system damping 
is achieved when auxiliary feedback is introduced in the control 
scheme. The stabilizing signal essentially modulates the 



reactive power in the system. 


1.4 Objftctlvftc and Scop« of This TKobIb 

This thesis is mainly concerned with the coordinated 
application of P.S.S. and S.V.C. for the improvementof dynamic 
stability of the power system. The main objectives of this 
thesis may be summarised as follows: 

1. To develop a suitable model of the power system in its 
simplest form consisting of a generator connected to an 
infinite bus through a transmission line and the S.V.C. 
located at the mid point of the transmission line. 

2. To design a P.S.S. which gives satifactory performance over 
a wide range of operation, assuming that only S.V.C. 
voltage controller is in the system and the S.V.C. 
stabilizer is not present. 

3. To design a S.V.C. stabilizer which gives satisfactory 
system performance over a wide range of operation, assuming 
that P.S.S. is not present in the system. 

4. To study the system performance when both the S.V.C. 
stabilizer and P.S.S. designed independently are operated 
simultaneously, and to study the robustness of the system. 

Robustness is one of the main requirement of the 


stabilizer. A stabilizer designed at an operating 
operate satisfactorily over a wide range operating 


point 
points . 


should 


1.5 R«vl*w of Lltorature 

Extensive research has been reported in the literature 
concerning the application of P.S.S. for the improvement of 



dynamic stability C1, 2, 3, 4, 53. A major problem with P.S.S. 
design is that P.S.S. designed at one operating point may not 
give satisfactory performance at other operating points. To 
overcome this problem adaptive control techniques have been 
proposed by many researchers. Irving II6D describe a Modal 
Reference Adaptive Control (MRAC) scheme. Ghosh et al. C73 
develop'^^self i^brting regulator for adaptive stabilization. A 
simple adaptive P.S.S. using gain scheduling approach has been 
developed by Brown Boveri C83. 

Bandyopadhyay and Prabhu C93 have given a novel gain 
scheduling approach to adaptive P.S.S. Madhu and Prabhu C103 
have adopted roubustness consideration in P.S.S. design. 

Application of S.V.C. for the improvement of dynamic 
stability is a recent development. Kinoshita C113 has reported 
the improvement of dynamic stability by auxiliary modulation of 
S.V.C. reactive power. Padiyar et al C1E3 have demonstrated that 
significant improvement in dynamic stability can be achieved by 
auxiliary stabilizing signal. Varma II133 has made a detailed 
study on the effects of S.V.C. stabilization on dynamic stability 
of a power system with a long transmission line. Ramar et al 
[1143 have shown how dynamic stability can be improved using a 
stabilizer in S.V.C. control circuit. 

Though many works have reported the application of P.S.S. 
alone or S.V.C. alone to improve system dynamic stability, very 

I 

few research reports are available on the coordinated application 
of P.S.S. and S.V.C. Hamouda et al C153 have discussed the 
technical advantage of coordinating S.V.C. and P.S.S. for damping 
inertial and torsional modes of steam turbine generators, \ 


locating S.V.C. at generator bus. C.H. Cheng et al C163 have 
discussed the application of P.S.S. and S.V.C. for a multimachine 
power system. 

1.6 Location of S. V. C. 

Location of S.V.C. in the power system for improvement of 
dynamic stability is an important problem to be carefully 
studied. 

Hamouda et al C153 have discussed the dynamic stability 
enhancement by locating the S.V.C. at generator bus. Kinoshita 
C113, Ramar C143 and Varma C133 have discussed the dynamic 
stability problem by locating S.V.C. at the mid point of the 
transmission line, as it will boost the midpoint voltage which in 
turn increases the power transfer capability of the transmission 
line as discussed in C133. 

Considering all these aspects, we have, in this thesis, 
located the S.V.C. at the midpoint of the transmission line. 

1.7 IdAntlf ylng Approprlat«r Signals 

In the design of P.S.S. or S.V.C. stabilizer proper 
selection of input signal for stabilization is important. 

In case of P.S.S., feedback of rotor speed deviation appears 
to be proper since the objective of P.S.S. is to provide adequate 
damping for electromechanical oscillations. Other signals such 
as terminal voltage, armature current, active and reactive power, 
field current, accelerating power, frequency etc. can also be 
used as P.S.S. input signals. Larsen and Swan 1153 have analysed 
the effectiveness of speed, power and frequency signals for 



P.S.S. Anwar C17D has examined the effectiveness of various 
P.S.S. input signals. Kundur C18D has reported the practical 
experience of Ontario Hydro, Canada, in P.S.S. use and has 
commented on the effectiveness of various P.S.S- input signals. 

Several research reports have come out examining the 
effectiveness of various input signals for S.V.C. stabilizers. 
Kinoshita CHI] has analysed stabilization using active power 
auxiliary signal for S.V.C. Padiyar et al C12II have demonstrated 
the efficacy of bus frequency signal for auxiliary controller. 
Varma C13I] has compared the effectiveness of various signals such 
as power, frequency, reactive power in improving the dynamic 
stability of system with a long transmission line. He has also 
introduced a new auxiliary signal concept called 'Computed 
Internal Frequency' (CIF) and has demonstrated its superiority 
over other signals. 

When both P.S.S. and S.V.C. are used in the system, various 
combination of signals can be made. Hamouda et al C15II have used 
modal speed as control signal for both P.S.S. and S.V.C. 
stabilizer. Research work has to be done to establish the 
effectiveness of various combinations of signals when both P.S.S. 
and S.V.C. are used simultaneously. 

1.8 SummAry of tho Work Roportod in thic Thosls 

A chapterwise summary of the work done in this thesis is 
given below: 

Chapter 2 deals with the development of system : model in 
state variable format for dynamic 'stability study. The state and 
output equations corresponding to P.S.S. using power and speed 



input signals are developed. Similarly output equations 
corresponding to S.V.C. stabilizer using power and frequency 
signals are developed. A numerical example is also presented. 

Chapter 3 deals with the design technique for design of 
P.S.S. and S.V.C. stabilizer, and methodology for coordination of 
PSS and SVC. An example problem is next considered and we 
present for it numerical results and robustness regions with 
P.S.S. alone, S.V.C. alone and then with both P.S.C. and S.V.C. 
in the system. Results are also presented with different 
combinations of signals. 

Chapter 4 contains conclusions and suggestions for future 


work. 



CHAPTER 2 

SYSTEM MODELLING 


2.1 Int.roduct.ion 

In the previous chapter the importance of P.S.S. and S.V.C. 
in damping electromechanical oscillations in power system is 
described. For the purpose of dynamic stability analysis, 
linearized model (State Space Model) of the system is used and 
the dynamic behaviour of the system is studied with the help of 
eigenvalue analysis. 

The power system in its simplest form consists of a 
synchronous generator connected to an infinite bus through a 
transmission line. In the literature, various state space models 
of the system having varying degrees of accuracy are available. 
The popular Hef f ron-Phi 1 1 i p model C193 has been used by many 
researchers. This model is simple, and offers a good physical 
understanding of the synchronous generator, but does not include 
damper windings. Therefore the results obtained with this model 
are pessimistic, since a major contribution to system damping 
comes fom damper windings C43. Padiyar and Ramshaw EI203 have 
developed a generalized synchronous machine model for large power 
systema dynamic studies. This model can take into account any 
number of windings in the d and q axes of the rotor. 

We for our study in this thesis, have considered Heffron 
Phillip model, as it is adequate for the purpose of dynamic 
stability study involving mainly the slow electromechanical and 
excitor mode behaviour. 



2. 2 Dev«l opmen't of Modol 

The state space model of the system is developed on the 
similar lines as in C193. 

The synchronous generator is connected to an infinite bus 
through a transmission line and S.V.C. is installed at the 
midpoint of the transmission line as shown in Fig. 2.1. 

The following assumptions are made for developing the state 
space model of the system given in Fig. 2. is 

1. The synchronous generator has no damper winding. 

2. Governor-turbine dynamics can be ignored. 

3. T model of the transmission line is sufficient. The line 
charging susceptance is clubb^sjed with S.V.C. susceptance. 

4. The machine stator and external network are in quasi-steady 
state. 

5. -Saturation in generator is neglected. 


2. 2. 1 Generator Model 

A third order generator model is considered. In the 
equation given below the standard notations are used. 

The dynamical equation corresponding to generator field 
circuit can be written as C21II 






’^d' ^d 


( 2 . 1 ) 


The dynamical equations corresponding to the generator swing 


equation can be written as 
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Figure 2.1 Single -machine Infinite -bus system with S.V.C. 
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In the equations given above 6 and Li are in radians and 
radians/sec. The rest of the quantities are in P.U. 


2t. 2. Z Excitation Systom 

We are interested in the study of dynamic stability 
problem caused due to thristor type excitation system. This is 
represented by IEEE Type IS CZZ3 model. Fig. 2.2 gives a block 
diagram representation of this excitation system. 

The dynamical equation representing the excitation system 
is 


PE 


FD 
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( 2 . 6 ) 


where is the generator terminal voltage. It can be 

represented in terms of Park's voltages as 


where 


2 2 2 

Vf = V j + V 
t d q 


V. = - X. 


(2.7) 


( 2 . 8 ) 
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2.2.3 Static VAR Ck^inpansator CSVO 

S.V.C. is normally a combination of reactors and 
capacitors, controlled by thristor circuits. Different types of 
S.V.C. configurations are used in practice. They are mainlys 

1. Fixed capacitor and thyristor controlled reactor (FC-TCR> 

2. Thyristor switched capacitor and thyristor controlled 
reactor (TSC— TCR). 

S.V.C. configurations mentioned above are shown in Fig. 2.3. 
The operating principle of S.V.C. may be explained in simple 
terms as follows. The voltage at the bus where S.V.C. is located 
is continuously measured. If the voltage measured is more than 
the specified voltage, capacitive susceptance is presented to the 
system. If the voltage measured is less than the sp e cified 
voltage inductive susceptance is presented to the system. This is 
achieved by varying the thristor firing delay. This varies the 
average current flowing through the reactor and consequently 
makes the reactor to effectively act as a variable reactance 
device. The general control scheme of S.V.C. is shown in Fig. 
2.4. 

Me have considered a simple S.V.C. voltage regulator in our 
modelling. Fig. 2.5 gives the block diagram representation of 
the S.V.C. voltage controller. 

The dynamical equation representing S.V.C. voltage 


controller is 
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Figure 2.3 (a) FC-TCR Figure 2.3(b) TSC-TCR 

Fc - Fixed capaciter 
Tc - Coupling transformer 
Th - Thyristor switch 
Re - Reactor 
Cl j C 2 ) C 3 Capacitors 





Figure 2,4 General control echeme of S.V.C. 





fig. Bloclr Diagr<tiii R<‘pr it .i on of GVC Voltage 
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2.2.4 Linearized Model 

The linearized state space model of the powersystem 

isobtained by linearizing eqns. (2.1), (2.2), (2.3), (2.6) and 

(2.10). The non— state variables AI., AV* , AT and AV are then 

d t e m 

eliminated using the following expression, Cderivation is given 
in Appendix AT 
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where parameters b^^,...,b 22 are defined in Appendix A. 

The resulting fifth order system model is described by the 
state equation 


px = Ax + b u. + b_ u_ (2.12) 

— — — "I — ^ ^ 

where 

X = [ A« a 5 AEpjj AB j (2.13) 

u^ = stabilizing signal from P.S.S. (2.14) 

Ug = stabilizing signal from S.V.C. stabilizer 

and the matrices A, b. and b„ are as given below: 

-1 -2 


(2.15) 







0 

^3 

14 

^5 ■ 





^21 

^22 

^23 

0 

^25 


A 

= 


0 

^32 

0 

0 

0 

(2.16) 





0 

%3 

^44 

«45 




- 


0 

^S3 

0 

«55 



0 

0 

0 

0 

] 



(2. 17) 

^2 = 

0 

0 

0 

0 bj 




(2.18) 

The parameters 


■ ■ ■ 5^5’ 

and b 

in the 

above 

equation are 

defined 

in Appendix 

A. 





2. 2. G Development of Output 

Equations 



The output equation 

of the 

system 

in general can be 

written 

as 









yi 

= 





(2.19) 

where represents 

the input signal 

to PSS, 

and 




^2 

= 

Cg X + 

dg Ug 



(2.20) 

where represents 

the input signal 

to the 

S.V.C. stabilizer. 

The elements 


the 

matr i ces 

c^ and 

Cg depend upon the actual 

variables y^ 

and yg 

chosen . 





2. 2. 5. 1 

Output 

Equation Corresponding lo PSS 

Input 

Signals 


For 

PSS we 

have considered 

generator 

active power and 


rotor speed as control signals. 



Power Signal! 


Generator active power is given by 


P = 3 (V, I . + 1 ) 

e d d q <1 


(2.E1 > 


In p.u., Pg = T^. Pg and are p.u. three phase power and 

torque, respectively. Hence we can write, 

AP = AT (2.22) 

e e 

The output equation corresponding to power signal input PBS is 
then, 

] ^ 


yi = = [‘■m ° •’41 ° ‘’42 


(2.23) 


Hence, 



I b.„ 0 

b. . 

0 

b.„ 1 


[40 

41 


^2 J 


(2.24) 


where parameters b^^, b^^ and are defined in Appendix A. 


42 


Spaed SI gnal i 

Here we use change in rotor speed (A<^) as PBS inpuit signal. 
Since speed is one of the state variables in our system model, 
the output equation is directly obtained as 

y, = ] 



X 


(2.25) 



<E.26) 


Hence c 


^ 0 1 0 0 0 j 


B. 2, S. Z Output, Equation Cor- r esponding to S. V. C. Stabilizer 
Input Signals 

For S.V.C. stabilizer we have considered mid-line power 
and S.V.C. bus angular frequency signals. 

Mld-llne Power Signals 

The mid— line power is given by 


P =3V.l. +V I 
m I md d mq q 


(2.S7) 


Linearizing the above expression and expressing sind AI^ in 

terms of state variables, we get 


^'"m - ^*60 ^^q *^61 * ^62 


(2.28) 


Thus we can write. 




= [ 


*>60 0 ■‘61 O “62 


] 


(2.29) 


The parameters and b^g are defined in Appendix A. 



SrVC Bus Ar»gular Fr wqucirtcy Slgrtsl 


S.V.C. bus angular frequency is obtained as 


d 

dt 




(2.30) 


Linearizing the above expression. 


V 




mdo 


(p AV > 
mq 


mo 


V 


rno 


(2.31) 


After manipulation (see Appendix A) we can write the final 
expression for Af^ as 




**80 ^23 ^81 ^82 *^83 | - ‘*2'^'^se 


(2.32) 




Hence , 


= [■ 


*^2 1*^80 **27 **81 **82 




(2.33) 


parameters bgQ» *^27' **81' **82' **83 **2 defined in Appendix 

A. 


2. 3 Nurnsrlcal ExamplA 

The power system data that we have considered in this thesis 
are as in Ref. C113, except the excitation system data. The 
system is as shown in Fig. 2.6. 



The generator represented in Fig. 2.6 is the equivalent 
generator representing all the generators of a generating 
station. The generator capacity is 5000 MVA. The power is 
transmitted over a distance of 200 km by a 500-KV two double 
circuit transmission lines <4 lines in parallel, in all) and SVC 
is installed at the mid— point of the transmission line. The 
power system data is given in Appendix B. 

We have chosen the infinite bus as the reference bus. The 
S.V.C. bus voltage is specified as 1.02 p.u. Assuming that the 
generator is operating at rated capacity at 0.9 power factor 
lagging, load flow is run to compute generator termj^inal 
conditions and the S.V.C. bus voltage angle. 

The A, and bg matrices of the system at the above 
operating point are 


■ -0.5283 

0 

-0.277 

0.156 

0.0296 ■ 

-106.60 

0 

-163.76 

0 

-2.465 

0 

1.0 

0 

0 

0 

-1212.14 

0 

- 9.64 

-100.0 

-408.98 

-59.07 

0 

14.48 

0 

-148.76 


b^ = O 0 O 5000 O 

bg = fo 0 O 0 1000 


C2.35) 

C2.36) 


The open loop eigenvalues corresponding to the above A 
matrix are 


-148.54 


-98.85 


-1.987 


-0.1116 ± j 12.79 




Fig. E.6 E 






2.x Concluding RomArkc 

In this chapter, we have developed the state-space model of 
a single machine infinite bus system with SVC located at the 
midpoint of the transmission line. The output equations 
corresponding to various input signals of the stabilizers are 
developed. A numerical example is also presented. 



CHAPTER 3 

DESIGN OF STABAILIZERS 


3.1 In-Lroduct^lon 

The state space model of the system developed in Chapter 2 
is used to design the P.S.S. and S.V.C. stabilizers. 

In this chapter P.S.S. and S.V.C. stabilizers are designed 
using the technique described in C233. We have adopted partial 
pole placement technique. The choice of assigned closed loop 
eigenvalues, design methodology and sector criterion are 
discussed in this chapter. Numerical results with different 
signals are presented. 

3.2 Problem Formulation 

The main objective of either P.S.S. or S.V.C. stabilizer is 
to provide adequate damping for electromechanical oscillations, 
de Mellow and Concordia C3T have provided on insight into the 
dynamic stability problem of synchronous generators, in terms of 
synchronizing and damping torques. Synchronizing torque is the 
component of the electrical torque in phase with rotor angle 
phasor and the damping torque is the component which is in phase 
with rotor velocity phasor. For the system to be dynamically 
stable, both these torques should be positive. Thus improving 
dynamic stability means providing an additional damping torque 
component, inphase with rotor velocity. This forms the basis for 
the design of P.S.S. or S.V.C. stabilizer as a phase lead network 
as shown in Figs. 3.1 and 3.2 respectively. 



Dynamic compensator 



Fi» 3 « 3*2 Structure of SVC Stabilizer 




Neglecting the filler and washout stages 


the P.S.S. and 


S.V.C. stabilizers have essentially second order transfer 
functions as given below: 


P. s. s. 


Fp(S) = 


(1 ^ STp^) <1 ^ STp3, 

(1 + STpg) (1 + STp^) 


or 


Fp(S) 




(3.1) 


S.V.C. Stablllzttri F (s) 
c 


(1 + ST ,) (1 + ST „) 
cl c3 

(1 + ST „) (1 + ST .) 

c2 c4 


or 


Fj.(S) = 


B s + d. s + B 
oc 1c 2c 


(3.2) 


3.3 I>esign Technique 

We make use of the output feedback pole placement algorithm 
of Munro and Hirbod C233, presented in Appendix C for P.S.S. and 
S.V.C. stabilizer design. The P.S.S. or S.V.C. stabilizer design 
uses linearized model of the power system since improvement of 
system response to small disturbances is of main interest. From 
the control engineering point of view, the role of P.S.S. or 
S.V.C. stabilizer can be viewed as shifting real part of 
eigenvalues corresponding to electromechanical mode to more 
suitable locations in the complex s-plane. 



Several techniques for design of stabilizers are available. 
But pole placement approach has the advantage that as in the 
classical root locus procedure, it gives the designer a feeling 
for system response. The algorithm of pole placement 
technique is simple. 

The algorithm discussed in ZZ33 provides a useful technique 
for the design of compensators for linear multivariable systems. 

If n be the order of the open loop system, q be the number 
of poles (eigenvalues) to be placed, m be the number of inputs 
and p be the number of outputs, the problem is to find an output 
feedback law of the form 

U(s) = - F(s) Y(s) (3.3) 

where F(s) is a proper, rational, polynomial matrix of order r, 
such that q poles of the closed loop system*., get placed at 
prespecified locations. If q < (n + r), the problem is referred 
to as partial pole placement problem. If q = (n + r), it is full 
pole placement problem. 

The order of the compensator r is related to q as 

q — max -C m,p > 
max i m , p J 

Here C □ means the next integer nearest to the quantity in 
brackets . 

In our case n = 5, r = 2 and max fm, p> = 1 since we are 

considering PSS and SVC separately. Hence the number of closed 
loop poles that can be assigned as per eqn. (3.4) is 3. During 
the design of stabilizers, we have observed that if we properly 



(3.4) 





prespecify three closed loop poles, the remaining unassigned 
poles will be at satisfactory locations in the left half of 
s-plane . 

In case of a single input P.S.S. or single input S.V.C. 
stabilizer, the transfer function matrix of eqn. (3.3) is a 
scalar function and has the form given in eqns. (3.1) and (3.2) 
respectively. 

3.4 Secbor Criterion, for Determining Robustness 

From the open loop eigenvalues presented in sec. 2.3 we 
observe that there is a pair of complex eigenvalues which 
correspond to the electromechanical mode of oscillation of the 
system. The real part of these eigenvalues is close to the 
imaginary axis in s-plane indicating that there is insufficient 
damping in the system. The main objective of the stabilizer is to 
shift these eigenvalues to a more stable location in the s— plane. 
The damping provided depends upon the closed loop locations of 
these eigenvalues. 

The region of robustness has been defined with respect to 
the sector shown in Fig. 3.3. The system is said to remain 
robust with respect to parameter variations if its eigenvalues 
remain within the sector. Note that the damping ratio of the 
dominant eigenvalues will be better than 0.0735. This damping is 
adequate for a powersystem. 

3. 5 Design and Coordination Methodology 

The methodology for the design of P.S.S. and S.V.C. 
stabilizer is essentially the same. At a design point three 




e Permissible Location of Close 


p.S.S. and three S.V.C. stablfi 1 izers are designed independently. 
Then a best pair of P.S.S. and S.V.C. stabilizer is chosen for 
coordination from robustness point of view. 

We present the design procedure adopted for P.S.S. design in 
detai 1 . 


3.6.1 P. S. £. D«rBlgn 

We have chosen a design point of 1 p.u. MVA loading at 0.9 
p.f. lagging. We have chosen this design point as this is in the 
normal rated operating range of synchronous generator and dynamic 
stability under heavy loading condition is of great concern. 
Furthermore, the system is expected to operate near this 
condition most of the time. At this design point a P.S.S. is 
designed assuming that S.V.C. stabilizer is not in the system, 
using the partial pole placement algorithm given in £ 231 . As 

already mentioned under Sec. 3.3, for an open loop system of 
fifth order only three closed loop poles can be placed at the 
desired locations using a second order stabilizer. The P.S.S. 
denominator is prespecified such that the stabilizer poles are at 
stable locations. The numerator of P.S.S. transfer function is 
determined by prespecifying three closed loop eigenvalues to be 
within the sector of Fig. 3.3. Two of the eigenvalues 
prespecified correspond to electromechanical mode of oscillation 
of the system. The locations of the unassigned closed loop 

eigenvalues are then checked and it is ensured that they lie 

within the sector of Fig. 3.3 

Similarly a few more P.Si$-are designed (we have actually 

designed two more PSS) at the same operating point by 



prespecifying three eigenvalues as in the first design explained 


above and choosing their locations to be different from one 
another but close to the corresponding location in the first 
design . 


3.5.2 S. V. C. Stabilizer Eieslgn 

After completing the P.S.S. design, several S.V.C. 
stabilizers (we have actually designed three stabilizers) are 
designed for the same design point and the same prespecified 
eigenvalue locations as for PSS design. The same procedure as 
discussed under PSS design (Sec. 3.5.1) is followed. 


3.5.3 Coordination of PS. S. and S.V.C. Stablllzar 

The design of P.S.S. and S.V.C. stabilizers is complete. 
The region of robustness of each P.S.S. designed is determined 
assuming that S.V.C. stabilizer is not present. Similarly the 
robustness region of each S.V.C. stabilizer is determined 
assuming that P.S.S. is not in the system. The procedure for 
determining the region of robustness is as follows. The operating 
point of the generator is changed in all the directions of the 
design point in P— Q plane (where P is the generator real power 
and Q is the generator reactive power) keeping the stabilizer 
parameters unchanged. The system parameters are determined for 
each operating point. The region of robustness of the stabilizer 
is that region where all the system closed loop eigenvalues 
remain within the sector defined in Fig. 3.3 

After completing the process of finding robustness of all 
stabilizers designed, it is observed that the region of 



robustness of P.S.S. is normally larger than the region of 
robustness of S.V.C. stabilizer. Now a trial and error method is 
adopted to coordinate P.S.S. and S.V.C. stabilizer. First a 
P.S.S. is chosen, correspondingly a S.V.C. stabilizer is chosen. 
Both of them are switched on to the system simultaneously. The 
region of robustness of the pair is determined. Similarly 
robustness region of all the other possible pairing of P.S.S. and 
S.V.C. (actually eight possible pairs in our case) are 
determined. After completing the process, it is observed that a 
further improvement in the region of robustness is obtained if 
the P.S.S. and the S.V.C. stabilizer paired have their individual 
region of robustness overlapped by largest amount. On the other 
hand if the P.S.S. and S.V.C. stabializer paired have their 
individual area of robustness quite different, the region of 
robustness of the pair will not be appreciable. It is normally 
less than the individual robustness of P.S.S. 

3.6 Choice of' Closed Loop Eigenvalues 

The open loop system has five eigenvalues. Out of these 
five two are in very stable locations.. The remaining three 
eigenvalues, two corresponding to the electromechanical mode 
(0.116 ± j12.79) and the third corresponding to exciter mode 
(-1.987), are near to the imaginary axis. So we have chosen 
these three eigenvalues for prespecifying. The real part of the 
eigenvalues corresponding to electromechanical mode are shifted 
near to -1.5 and their imaginary parts are changed depending upon 
the requirement. The eigenvalue corresponding to the exciter 
mode is shifted near to -2.0. The exact locations of the 



prespecified eigenvalues is determined by an computer aided 
trial and error method keeping in mind that the unassigned poles 
should also lie within the sector defined in Fig. 3.3 and 
stabilizer gains should not be too high. 

3. 7 Numerical Results 

For all the numerical results presented in this section, the 
design point is that corresponding to 0.9 p.f. lagging on 1 p.u. 
MVA semicircle. The open loop ei g env a lues of the system 
corresponding to this operating point are given in section 2.3. 

3.7.1 Results with P. S. S. 

Here we present the numerical results obtained with PSS in 
the absence of SVC stabilizer. Results are presented for power 
input and speed input PSS. 

3.7.1. 1. PSS with OeiMrator Active Power Deviation AP as Input 
In this case, 

y-i = ^ 

= Ko ° *>41 ° ‘■42 ] <3.5) 

where the parameters are defined in Appendix A. 

Table 3.1 presents the numerical results pertaining to three 
PSS designed at the same design point. The result presented 
shows the closed loop eigenvalues assigned, transfer function and 
the range of robustness of each PSS on 1 p.u. MVA semicircle. 



).7. 1>2 PSS wlt-h Angular Valocit,y Daviai/lon A<-i Ab Input 
In this case, ’ 

= A<^ = X 

= [0 1 O O olx (3.6) 


Table 3.S presents the closed loop eigenvalues assigned, 
transfer function and range of robustness on 1 p.u. MVA 
semicircle of each PSS designed. 


7. 2. Results with SVC Stabilizer 

Here we present the numerical results obtained with SVC 
stabilizer in the absence of PSS. Results are presented for 
power input and frequency input SVC stabilizer. 

3. 7. 2.1 SVC Stabilizer with Midline Active Power AP^ as Input 
In this case. 


= AP 


m 

Cg X 



lb. „ 

0 


0 

b, _ 1 


[60 


61 


62 J 


(3.7) 


where the parameters are defined in Appendix A. 

Table 3.3 presents the numerical results pertaining to three 
SVC stabilizers designed at the same design point as for PSS 
design. The result presented shows the closed loop eigenvalues 
assigned, the transfer function and range of robustness of each 
SVC stabilizer on the 1 p.u. MVA semicircle. 



AS 


3.7"»££«^ SVC St/AbilizAr wl'th SVC Bus FrAquAncy DAviAbion A<4> 
Input 


In this case, 

yg = A, 


“ '^2 S * <*2 “j 


= [ 


^80 *^27 *^81 *^82 


"83 




d„ AV 
2 se 

<3. 8) 


where the parameters are defined in Appendix A. 

Table 3.4 presents the numerical results pertaining to three 
SVC stabilizers designed at the same design point. The result 
presented shows the closed loop eigenvalues assigned, the 
transfer function and range of robustness of each SVC stabilizer 
on the 1 p.u. MVA semicircle. 


3.7.3 Numerical Results with both PSS and SVC Stabilizer 
Coordi ruaited 

PSS designed under section 3.7. 'I and SVC stabilizer 
designed under section 3.7.2 are coordinated according to the 
methodology outlined in section 3.5. Results are presented when 
PSS and SVS stabilizers and simultaneously operated as given 
below: 

(i) AP input PSS and AP input SVC stabilizer 

(ii) Aa input PSS and Aw^ input SVC stabilizer 

(iii) AP^ input PSS and A6i>^ input SVC stabilizer 
A6> input PSS and AP^ input SVC stabilizer 


(iv) 



3. 7. 3.1 AP Input PSS and AP Input SVC Stabilizer 
^ m 

Table 3.5 presents the results for the coordinated 

application of AP input PSS and AP input SVC stabilizer. The 

e m 

result presented shows the closed loop eigenvalues, transfer 
function of PSS and SVC stabilizer chosen for coordination and 
the range of robustness of 1 p.u. MVA semicircle. 

Figure 3.4 shows the region of robustness for PSS alone, SVC 
stabilizer alone and for coordinated PSS and SVC stabilizer 
case. 


Figure 3.5 shows the initial condition response of A6 in all 
the three cases. 


Table 3.5 presents the range of robustness 

on 1 

p. u. 

MVA 

seicircle of all 

other pairing 

possible 

in this 

case 

which 

are 

not selected for 

coordination. 

and Fig. 

3.6 shows 

the 

region 

of 


robustness for all such cases. 

3. 7. 3. 2 JUji Input PSS and Aoi^ Input SVC Stablllzar- 

Table 3.7 presents the results for coordinated 

application of A6> input PSS and A6i^ input SVC stabilizer. 

Figure 3.7 shows the region of robustness on P-Q plane. 
Figure 3.6 shows the initial condition response of A<5 . 

3.7.3. 3 AP^ Input PSS and Input SVC Stablllzar 

Table 3.9 presents the results for coordinated 

application of AP input PSS and A<.> input SVC stabilizer. 

© s 

Figure 3.10 shows the region of robustness on P-Q plane. 
Figure 3.11 shows the initial condition response of A6 . 
Table 3.10 presents the range of robustnessa on 1 p.u. 



MVA semicircle of other possible pairing in this case which are 
not selected for coordination and Fig. 3.12 shows the region of 
robustness for all such cases. 

3. 7. 3. -4 A<.> Input PSS and AP Input SVC S^ablllzar 

Mi 

Table 3.11 presents the results for coordinated Ato input 

PSS and AP input SVC stabilizer, 
m 

Figure 3.13 shows the region of robustness on P-Q 

plane . 

Figure 3.14 shows the initial condition response of A<5 . 
Table 3.12 presents the range of robustness on 1 p.u. 
MVA semicircle which are not selected for coordination and Fig. 
3.15 shows the regioan of robustness for all such cases. 

3.7.4 Discussion of Results 

After a careful study of design and coordination results, 

we are able frame »the following set of rules to coordinate 

A 

P.S.S. and S.V.C. stabilizer in a system. 

1. Choose a design point near the rated range of the generator 
and design several P.S.S. and S.V.C. stabilizers 
independently, by pole assignment. 

2. Determine the region of robustness of each P.S.S. and S.V.C. 
stabilizer independently. 

3. Select that P.S.S. and S.V.C. stabializer for pairing whose 
individual regions of robustness are large and overlap by 
the maximum amount. The resulting pair will further improve 
the system performance significantly. 



Tabl^ 3.1 Numerical Rasulla for Inpurt PSS 


SI. No. 

Closed loop PSS transfer function 

eigenvalues 

assigned 

Range of 
robustness on 

1 p.u. MVA 
semicircle 
(P,Q) 

1. 

-1.54±j.13.Z F (s) = 0.098s^+3.40s+Z3.6 

(0.83,0.56) 


® s^ + Z4s + 80 

To 


-2.5 

(0.816,-0.57) 

Z. 

- c;* ^-5 « c / . 0.114s^+3.Zs+16.6 

-1.5±j13.0 F (s) = R 

® s^ + 24s + 80 

(0.842,0.54) 


To 


CVJ 

1 

CVJ 

1 

(0.83, -0.56) 

3. 

-1.5±J12.a F (s) . 

(0.76, 0.65) 


® s^ + 24s + 80 

To 


-2.0 

(0.876,-0.78) 





Tabla 3. 2 I Numcirlcal Rcrsul'ts for A(i> Inpixt PSS 


SI. No. 

Closed loop 
eigenvalues 
assigned 

PSS transfer function 

Range of 
robustness on 
1 p.u. MVA 
semicircle 
(P,Q) 

1. 

-2.5±j13.2 

-2.0 

P , , -0.138s^-0.82s-0.51 

F (s> = 

® s^ + 30s + 200 

(0.852, 0.523) 
To 

(0.948,-0.324) 


-2±j13.0 

-2.0 

^ , -0.113s^-0.426s-0.22 

C. # «- 1 — 

(0.82,0.572) 

To 

(0.96, -0.28) 

CL m 

^ ^ O 

s^ + 30s + 200 


-1.5±j13.6 

- 2.0 


F (s) = 
s 


-0 . 092s -O . 8 1 s- 1 . 07 
+ 30s + 200 


(0.858,0.513) 

To 

(0.902, -0.43) 




Tabl« 3.3 I Numerical Recults for AP Input SVC SJtablllzer 


>l.No. Closed loop 

eigenvalues 
ass i gned 

S.V.C. stabilizer 
transfer function 

Range of 
robustness on 

1 p.u. MVA 
seiiiicircle 
(P,Q) 

1. -1.5±j13.0 

-2.2 

^ , -0.2535^-4.45-10. 51 

F (s) - 

s"^ + 24s + 80 

(0.840,0.542) 

To 

(0.952,-0.30) 

1. -1.54±j13.-1 

-2.5 

_ , , -0.213s^- 5. 136S-15.41 

(0.864, 0.50) 

To 

(0.946,-0.324) 

^ ^ p 

4 24s + 80 

3. -1.5±j12.8 

-2.0 

r. , , -0.271s^-4. 0865-7. 3 

F (s) - 

s^ + 245 + 80 

(0.862, 0.506) 
To 

(0.92, -0.40) 




TAbl« 3. 4 I NumArlcAl ReBult.s f‘or Au> Input, SVC SLablllzer 

B 


Sl.No. Closed loop 

e i gen val ues 
assigned 

S.v.c. stabilizer 
transfer function 

Range of 
robustness on 
1 p.u. MVA 
semi circle 
(P,Q) 

1. -1.56±j 10.93 

-2.2 

-0.0315-2. 75-6.73 

(0.861,0.51) 

To 

(0.944, -0.33 

5 

5 + 24s + 80 

• 

-0.027s^-2.67s-5.31 


2. -1.5±j10.89 F (s) = 

5 

(0.88,0.475) 

c 

5*^ + 245 + 80 

To 

-2.0 


(0.92, -0.39) 


3. -1.54±j13.2 

-E.5 


F (s) 
c 


-0.046s -2.85 - S.-q 
5^ + 245 + 80 


(0.834,0.55) 

To 

(0.962, 0.273) 
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a/ 

Tables 3 *^ Mum»rlcal Rttsulbs for Coordlnabod AP 

o 


and AP Incut SVC Stabilizer 

m 


Transfer functions of PSS 
and SVC stabilizer chosen 
for coordination 


Location of 
closed loop 
eigenvalues 


PSS 


Fs<s) = 


0.08s +2.89S+6.3 


— 2 

s + E4s 


80 


SVC stabilizer 


F^(s) = 


-0.271s -4.08S-7.3 
S+24S+80 


-140.2 


- 78.3 

- 3.14±j12.4 

- 22.59 

- 2.0 

- 5.22 

- 20.23 

- 3.99 


I nput PSS 


Range of 
robustness on 
1 p.u. liVA 
semi ci r cle 
(P,Q) 


<0.7, 0.71) 

To 

(0.86, -0.51) 



ST 


IV- 



nd d 


Fig. 3.4 Region of Robustnes for AP Input PSS and AP Input 
SVC Stabilizer ® 




PSS, AP Input SVC Stabilizer 
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T*ble 3.6 Numerical Results for Other Pairs of input P.S.S. 

and Af^ input S.V.C. Stabilizer not 
Selected for Coordination 


SI .No of P.S.S 
chosen from 
Table 3 . 1 

SI .No of S. V. C. 
Stabilizer chosen 
from Table 3.3 

Range of 
robust ness 
on 1 p . u . 

MVA semicircle 
(P , Q) 

Sub-F i g . No . 
under 

Fig. No 3.6 

1 

1 


(0.83 , 0.56) 
to 

(0.8£ ,-0.57) 

Fig. ( a } 

1 

£ 


(0.84 , 0.54) 
t o 

(0.83 , -0.56) 

Fig. (b) 

1 

3 


(0.76 , 0.65) 
t o 

(0.86 , -0.48) 

Fig. (c) 

S 

1 


(0.84 , 0.54) 
to 

(0.98 , -0.££) 

Fig. (d) 

£ 

£ 


(0.84 , 0.54) 
to 

(0.965 , -0.£6) 

Fig. ( e ) 

2 

3 


(0.84 , 0.54) 
to 

(0.91 , -0.4£) 

Fig. (f ) 

3 

1 


(0.80 , 0.6) 
to 

(0.91 , -0.4£) 

Fig . ( g ) 

3 

£ 


(0.80 , 0.6) 
to 

(0.89 , -0.46) 

Fig . ( h) 














Table 3.7 Numerical Results for Coordinated 
Aci)^ Input SVC Stabilizer 


Transfer functions of PSS 
and SVC stabilizer chosen 
for coordination 


Location of 
closed loop 
eigenvalues 


PSS 

Fg(s> 


-O . 1 13s^-0 . 426S-0 . 22 

2 

s + 30s + 200 


-249.5 


- 95.58 

- 22.69 

- 3.18±j11.5 

- 5.26 


SVC stabilizer 


F^(s> 


-Q.046S -2.as-84 
s^ + 24s + 80 


4.82 

2.48 


Input PSS and 


Range of 
robustness on 
1 p.u. MVA 
semicircle 
(P,Q) 


(0.782, 0.62) 
To 

(0.942,-0.33) 




Region of Robustness for As. Input PSS and As, 

Input SVC Stabilizer 
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Table 3.6 Numerical Results for Other Fairs of AcJ input P . S . S . 
and input S.V.C. Stabilizer not 

Selected for Coordination 


SI . No of P . S. S 
chosen from 
Table 3,2 

SI .No of S, 
Stabilizer 
from Table 

V. C. 
chosen 

3 . 4 

Range of Sub-Fig. No. 

robustness under 

on 1 p. u . Fig . No 3 . 9 

HVA semicircle 
(P , Q) 

1 

1 



(0.65 , 0.54) 
to 

(0.96 ,-0.28) 

Fig . ( a ) 

1 

2 



(0.85 , 0.54) 
t 0 

(0.95 -0.30) 

Fig . (b) 

1 

3 



(0.65 , 0.5) 
to 

(0.95 , -0.3) 

Fig. ( c ) 

£ 

1 



(0.846 . 0.53) 
to 

(0.96 , -0.28) 

Fig. (d) 

2 

2 



(0.858 , 0.513) 
to 

(0.96 , -0.28) 

Fig. (e) 

3 

1 



(0.858 , 0.513) 
to 

(0.946 , -0.324) 

Fig. ( f ) 

3 

£ 



(0.858 , 0.513) 
to 

(0.946 , -0.324) 

Fig . { g ) 

3 

3 



(0.858 , 0.513) 
to 

(0.98 , -0.19) 

Fig. (h) 


11 224 . a 


Fig . 



Q pu Plg<<i) Q 









Fig(g) Q PU Flg(h) 
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Table 3^ 


Numerical Results for Coordinated AP^ Input PSS 
and Aoi Input SVC Stabilizer 

fi 


Transfer functions of PSS 
and SVC stabilizer chosen 
for coordination 


Location of 
closed loop 
eigenvalues 


Range of 
robustness on 
1 p.u. MVA 
semi circle 
(P,Q> 


PSS 


F^(s) 

s 


0.08s^+2.S9s+6.3 

2 

s + 24s + 80 


SVC stabilizer 

^ , -0.027s^-2.67s-5.31 

F (s) = 2 

s'^ + 24s + 80 


-176.3 


- 97.3 

- 21.8 


(0.744,0.67) 

To 


- 3.76±j11.46 (0.84,-0.54) 

- 6.46 


- 4.2 


- 2.49 


3.98 






Time in Sec 



Initial Condition Response of A«5 ; Input 
PSS, Aw Input SVC Stabilizer 












•'‘S''” 0 pu 
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Table 3J-1 : Numerical Results for Coordinated Ato Input PSS and 

AP^ Input SVC Stabilizer 


Transfer functions of PSS 
and SVC stabilizer chosen 
for coordination 

Location of 
closed loop 
eigenvalues 

Range of 
robustness on 

1 p.u. liVA 
semicircle 
(P, Q) 

PSS 

-139.87 


^ , -0.09Es^-0.ais-1.07 

F (s) = 2 

® s + 30s + 200 

- 82.83 

- 3.48±j12.42 

(0.824, 0.56) 

To 

(0.898, -0.44) 

SVS stabializer 

-12.99 


^ , , 0.271s^-4.0S6s-7.23 

F (5) = 2 

^ s + 24s + SO 

- 2.01 
- 4.75 
-19.36 





Pig* 3.13 Region of Robustness for Input PSS and 
PSS and Af> Input SVC Stabilizer 
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Table 3.18 Nuaerical Result* for Other Pairs of Ago input P.8.8, 
and input S.V.C. Stabilizer not 

PT) 

Selected for Coordination 


81. No of P.8. 8 
chosen froa 
Table 3.8 

81 .No of S.V.C. 
Stabilizer chosen 
from Table 3.3 

Range of 
robustness 
on 1 p . u . 

MVA semicircle 
(P , Q) 

Sub-Fig. No. 
under 

Fig. No 3.15 

1 

1 

(0.85 , 0.58) 
to 

(0.96 ,-0.878) 

Fig. (a) 

t 

8 

(0.88 , 0.47) 
to 

(0.968 , -0.878) 

Fig. (b) 

1 

3 

(0.88 , 0.47) 
to 

(0.968 , -0.878) 

Fig. (c) 

8 

1 

(0.88 , 0.578) 
to 

(0.96 , -0.88) 

Fig. (d) 

8 

* 8 

(0.88 , 0.57) 
to 

(0.98 , -0.88) 

Fig. (e) 

8 

3 

(0.88 , 0.57) 
to 

(0.98 , -0.8) 

Fig. (f ) 

3 

1 

(0.86 , 0.51) 
to 

(0.98 , -0.41) 

Fig. (g) 

3 

8 

(0.86 , 0.5) 
to 

(0.90 , -0.43) 

Fig. (h) 


ct#d for coord 


















3.8 A F®w Rttmarks on Doslgn 


While 

determining 

the 

region 

of robustness we 

have 

maintained 

the S.V.C. 

bus 

voltage 

constant at 1.02 

p.u. 


irrespective of the power being transmitted in the line. This 
can call for a higher capacity of S.V.C. to be installed. But in 
practice, S.V.C. voltage controller reference setting is made to 
change depending upon the power flowing through the S.V.C. bus. 
This is to optimize the rating of S.V.C. to be installed. To 
simplify our design procedure we have maintained a constant 
voltage as mentioned above. The validity of our results have to 
be reestablished for such practical situations where the S.V.C. 
reference voltage is changed with respect to power flow. 

The second remark is on the power factor of the generator, 
when the S.V.C. is operating at the midpoint of the transmission 
line with the instruction to maintain certain level of voltage, 
the generator power factor is mainly decided by the power demand 
and S.V.C. rating. This is quite a encouragiuj, factor to the 
results we have presented as it really implies that in practice 
the stabializer that we have designed will not encounter most of 
the operating points in the robustness region at all. Our 
intention was to highlight the fact that under any operating 
condition of the generator the scheme that we have proposed can 
work well. 

Finally a few words about the region of robustness. We have 
tested the robustness range of the stabilizer by changing the MVA 
output of the generator from 1.E p.u. MVA to 0.2 p.u. MVA in 
steps of 0.1 p.u. MVA. 



Bui inpractice, the stabilization may not be required when 
the generator is working at low loads. Dynamic stabi^lity 
problem is usually encountered near full load operation. At low 
loads even if the low frequency oscillation exists, it is not as 
serious as at high loads, and it may not lead to the tripping of 
the t^ransmission lines etc. as the transmission lines are 
operating far below their rated transmission capacity. 

3.0 Concluding Romarkc 

In this chapter we have designed PSS and SVC stabilizer for 
a power system. Stabilizers are designed considering power and 
angular frequency as input signals. Then we have coordinated both 
the stabilizers by computer aided trial and error method. At the 
end of this process we have come up with a set of empirical rules 
to coordinate P.S.S. and S.V.C. stabailizer to give better 
system performance. Robustness criterion has been chosen to 
grade the performance of a stabilizer. From the results 
presented in this chapter, it is also clear that if one of the 
stabilizers fail to operate, the system is robust enough in the 
individual robustness range of the stabilizer which is working. 



CHAPTER 4- 
CONCLUSIONS 

In this thesis we have shown how P.S.S. and S.V.C. 
stabailizer designed independently can be coordinated to obtain a 
further improvement in system performance from dynamic stability 
point of view. 

We have presented a new model of the power system for the 
dynamic stability analysis of a generator connected to an 
infinite bus through a transmission line with S.V.C. at the 
midpoint of the transmission line. The model proposed is simple 
and straight forward. 

P.S.S. and S.V.C. stabailizers are designed independently 
using pole placement technique. We have then coordinated the 
P.S.S. and S.V.C. stabilizer by computer aided trial and error 
method. In all the design and coordination processes robustness 
was our main concern. At the end of the coordination process we 
have come up with a set of rules to coordinate P.S.S. and 5^.V.C. 
stabilizer in power system. 

The coordination process is done considering power and 
speed signal for P.S.S. and power and angular frequency signal 
for S.V.C. stabilizer. The results presented clearly indicate 
that in all the combinations of P.S.S. and S.V.C. that we have 
chosen for coordination, if one of the stabilizer fails to 
operate, a minimum region of robustness equal to the individual 
robustness of the stabilizer that is working is always ensured. 



SUGGESTIONS FOR FURTHER WORK 


1. We have not fixed the rating of S.V.C. in this thesis. It 

has been assumed that the S.V.C. has enough capacity to 

maintain the voltage that we specified at S.V.C. bus. But in 
practice there is a limit for S.V.C. operation both on 

inductive and capacitive range, fixed by its rating. 
Studying the system performance fixing SVC rating is an 
interesting area to work. 

2. We have considered a simple third order generator model, 

neglecting damper circuit. The results obtained in this 

thesis may be tested on a higher order generator model 
taking into account damper windings. 

We have considered single input stabilizer. Work 
carried out considering two input stabilizers. 


3. 


can be 
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APPENDIX A 

STATE SPACE MODEL OF THE SYSTEM 


A. 1 Etevelopment of the Model 

The s/stem model is developed on lines similar to those of 

C19D. 

The power system consists of a synchronous generator 
connected to an infinite bus through a transmissioan line. 
Static VAR compensator (SVC) is installed at the midpoint of the 
transmission line (See. Fig. <2.1>). 

The following assumptions are made for developing a state 
space model of the system. 

1. The synchronous machine has no damper windings. 

2. Governor and turbine dynamics can be ignored. 

3. The machine stator and external network are in quasi— steady 
state . 

4. Magnetic saturation in generator is neglected. 


A. 1 . 1 Generator Model 

A third order generator model is considered. The equation 
and notations are standardr see for example, C21I]. The equations 
are given below: 




^d' ^d 


Pw 


(T 


m 


- T ) 
e 



(« — <0 ) 
o 


(A.1) 


M 


(A. 2) 



P6 


(A. 3) 


01 - Oi 

o 

where 

M = ?H_ 

6> 

O 

= 3<V.I. + V 

0 da 4 4 

V . = -X I 

d q q ^ 

V = x’ I . + E* CA.4) 

q q d q 

All quantities are in p.u. except 6 and oi, which are in radian 
and rad/sec respectively. 


A. 1.2 Excitation System Model 

A single time constant excitation system, IEEE Type IS 
C223 is considered. Fig. (2.2) gives a block diagram 
representation of the system. 

The dynamical equation representing the excitation system is 


PE 




K. 


FD 


<^ref - ^t 




(A. 5) 


where can be represented in terms of Park's voltages as 

Ir 


V* = (V. 


+ v2,1/2 


where 


= — X I 

q q 

“ X . I . + E 

d d q 



A. 1 . 3 SVC Model 


A single time constant SVC voltage controller is 
considered. The block diagram representation of SVC voltage 
controller is given in Fig. (2.3). 

The dynamical equation representing SVC voltage controller 
is 


PB = 



+ 


— (U 

Tg ^'^mref 


m 


V ) 
se 


(A. 6 ) 


A. 1.4 Linearized Model 

The linearized state space model is derived by linearizing 

eqns. (A.1), (A. 2), (A. 3), (A. 5) and (A. 6 ). The non state 

variables AI AT , AV* and AV are represented in terms of state 
d e t m 

variables as follows. 

A. 1.4. 1 Expressions for AI . and A1 

d q 

To obtain expressions for AI . and AI in terms of state 

d q 

variables the system in Fig. (2.1) is reduced to the form shown 
in Fig. (A.1) by star-delata transformation. The impedance then 
reflecting between the infinite bus and the ground can be 
neglected. 

For the system in Fig. A.1, we can write the following 
equations as in C21II, Sec. 5.4, 

V. (1-»A .) + V X = - V sin ( 6 - 0 ) + R I + x I 
d 1 4 ^ ^ edeq 


(A. 7) 







-V + V (14X .) = V cos(6-a) - x„I. +R„I„ 
dZq1 «o edeq 


(A. 8) 


where ot and & are as defined in C21D, Sec. 5.4, 


«10 = « 


R = 2R - 2XB 
e 



O 

= (X 

2 2 
- x^ = R B 


^10 

2 2 2 
= R + (X 1-) 


X 

(R.^R +X.,,X ) 1 

_ 10 e 10 e ^ _ 


^1 


^2 ' 2 

10 

Linearizing equation 

(A. 7) and (A. 8) we < 

A1 . and Al in terms 
d q 

of state variables as 


’ ^"d ' 


**11 **12 **13 


Al 

q 


**21 **22 **23 


R._X - X._R 
10 e 10 e 
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(A. 9) 
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- 

AB 


(A. 10) 


where 


**11 


b- B - b_(1-B X) 

3 O £: O 


‘l **5 " 


** 2**4 


(A. 11) 


'12 


= - [- 


b,.V cos(<5„ - a) - b„ V - a) 

5 00 o g oo o 

**1 **5 “ **2 **4 
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(A. 12) 



(A. 13) 


'13 


L “i 


“5 ~ *^4 


'21 


-[- 


b. B - b^(1 - B X) 
4 o 1 o 

bg b^ - b^ b^ 


22 


= -[- 


b- V cos(6 —a) — b.V sin(6 — a) 
4 06 o 1 tut o 


“2 


bi b^ 


23 


f ~ 1 

•- b b. - b. b- J 
2 4 1 5 


In the above equations 


b. = B X . R - 2R + 2B RX 

1 o d o 

b„ = -X + B X X - R^ B + X^B - 2X 

2 q o q o o 


= -[ 


I X X + I . X .R + E R + 21 . RX 
qo q do d qo do 


I R^ + I 
qo qo 


*>4 = 


[=<4 - *d 


.X B + R^B - X^ B + 2X 
do o o 


] 


= [ 


X RB - 2R + 2 RX B 
q o o 


= -[> 


b^ = - |X RI - x’ X I - X e' + R^Irf„ - X^I. + 2RX 
6 I q qo d do qo do do 


<A.14) 


(A. 15) 


(A. 16) 




So] 


(A. 17) 



A. 1 . 4. a 


Expressions for AT 


where 


'*'0 = 3 + V I > p.u. 

e d d q q 


V . = - X I 

d q q 


V„ = X . I . + E 
q d d q 


Therefore T 


= " K 


^ I - (X - X .) 
q q q d d q 


(A. 18) 


Linearizing equation (A. 18) and then substituting for AI^ and 
from equation (A. 10) we get 


AT 


= m[ 


'‘40 + *>43 


(A. 19) 


where 


*>40 = ^[•‘zi'So - Xql 


X. - I X ) + I 

do d do 11 qo d qo q qo 


b.„ = 3^b^^ (E__ - X_ I^_ + X4 I^_) + b^„ (I__ X^ - I__ X 


'4E 


2Z qo q do d do iZ qo d qo d 


:>] 


*-43 = 3 


[•■hs'So - 


XI. + X . I . ) + b.., (I X . - I X . 
q do d do 13 qo d qo d 


-] 


(A. 20) 


A. 1.4. 3 Expression for AV^ 

For generator terminal voltage we can write 



(A. 21) 


vf = 

t d q 

Linearizing above equation and substituting for AI . and 

d 

Al^from equation (A. 10) we can write 


AV = K, AE + K- A«5 + K„ AB 

t o q 3 9 


(A. 22) 


where 


•<6 - 


do 

'to 


^q **21 V 


qo 


to 


^'d **11 


qo 


to 


s = 


do 

^ ""22 


go Y K 


K, 


do 


go y’ h 

'^to ** 


(A. 22) 


A. 1 . 4. 4 Expression for AV^ 

For SVC bus voltage V we can write 

m 


2 2 

V = V^ 
m mq 


md 


where 


(A. 23) 


V.= V. - I. R-I X 
md d d q 

V = V - I R + I .X 
mq q q d 


Linearizing equation (A. 23) and subsdtituting for ^nd 

from (A. 10) we get 


AV = + l>=-7 A<5 + b__ AB 

m 56 q 57 58 


(A. 24) 



The final state equations are obtained after substituting 
^"^e' ^'^^t terms of state variables, in eqns. 

(A.1), (A. 2), (A. 3), (A. 5) and (A. 6). The resulting fifth order 
system state equation can be written as 


X = A >1 + u-l + bg Ug 


where 


« =[K 


^ Aci AS AEpjj AB 


(A. 25) 


(A. 26) 


and u^ = ^'*^ss' stabilizing signal from PSS 

Ug = stabilizing signal from SVC stabilizer 


A^ , b^ and bg matrices are given as 


(A. 27) 
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where 
^11 ' 


C< X . - X . ) b., - 13 
d d 11 


do 



13 




do 


'*d - "d> *>13 


14 


*15 


do 


do 


= 3 (E__ - X_I,_ + X, _!,_)+ b^^(X^ I__-X_I__)+ I 


-r A. Xj I-r . X “ A X # -r X 

21 qo q do do do 11 d qo q qo qoj 




^ •T’A. X. /-rc iAjX —A X 

22 qo q do do do 12 d qo q qo 


^25 " ^ [**'’'* - ^ 


\ C. ~ A X j 

23 qo q do 


X* Ij ) + b^„(x' I -XI )1 
do do 13 d qo q qo J 


K, 


”22 

^2 

^1 = 
^43 


M 


^6 »<A 


S ‘^A 


”44 

^45 

Si = 
S3 = 


1 

^A 

S ‘<A 


S S6 
S 

•<8 S7 

S 


Ss = 


' ^ S Sa 


B 



CA.31) 



A. 1 . 5 Output Equations 

The output equation of the system can be written as 

yi = X 

where represents input to PSS, and 

Vs = "2 >1 

where y^ represents input to SVC stabilizer. 

The elements of matrices c^ and c^ depend upon the 
particular variables chosen as y^ and y^* 

A .1 . 5. 1 Output Equations Corresporwllng to PSS SI gruels 

For PSS we have considered rotor speed and generator 
accelerating power as control signals. 


S^teed Signal 

Chanage in rotor speed (A(i>) is utilized as control signal. 
Since speed is one of the state variables in our system model, 
the output equation is directly obtained as 


y^ = LAo>l = 1^ 0 
Thus, '•2 ~ f ^ 


10 0 


O J X (A. 32) 
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Power Signal 


Generator power output is taken as control signal. When 
expressed in p.u. AP^ = H^f'ce we can write the output 
equation corresponding to power signal as 


= [ ^% ] = [*> 


40 ^ ‘‘aI ° ^4Z\ ^ 


Thus, 


(A. 34) 



[b.. 

0 

b,^ 

0 

b . _ 1 


L 40 
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where b. 


40 “42 

(A. 19) and (A. 20)). 


b,., are as given in the expression for AT , (eqn, 


A. 1 . S. 2 Output E)quatlons Corresponding to SVC STtablllzer Signals 
For SVC stabilizer we have considered midline power and 
SVC bus angular frequency signals. 


Midline Power Signal: 

The midline power is given by 

^ 1^ I, ] < 

Linearizing the above expression and manipulating we get 


AP = b,^„ AE + b. . A«5 + b._ AB (A. 37) 

m AO q A1 A2 
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51 
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22 

mqo 


54 

qo 

62 = 

^3 

V . 
mdo 

+ 

^52 

^do 

+ 

*^23 

V 

mqo 

4“ 

S5 

I 
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In the above equation 

•>51 = [ '->‘q - *> ■>21 - « *>11 ] 

•>52 “ [ <-*, - >•' ^>22 - « ^>12 ] 

■>53 = [ - >" ^>23 - ^ ^>13 ] 

■>54 = [ <-'■■>21 ♦ <■< ^ >'d'^>1l"'' ] 

■>55 = [ '-’•■>22 ^ - ■'d>^‘l2 ] 

■>56 = [ <-'■■>23 * '* * ’'^'■>13 ] 


(A. 40) 



SVC Bus Angular Frequency Signal s 

For SVC bus angular frequency we can write 


<>> 


d 

dt 


r * ^md "1 

[ — I 

^ mq J 


Linearizing the above expression. 
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Linearizing (A. 42), 


AV 


m 


= [-^ ] - 

^ mo '■ mo ^ 


Also we can write 


V ^ = V . - I .R - I^X 

md d d q 


V = V - R + I^X 
mq q 4 “ 


(A. 41) 
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In terms of incremental quantities, we get 



(A. 44) 
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Substituting eqns (A. 43) appropriately in eqn. (A. 41) we get the 
final expression for A6>s as 


] = [b 
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Also, “ ^ 11 ' **36 “ ^13' **37 ^14' ®38 ^ 
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APPENDIX B 

POWER SYSTEM DATA 


Most of the data are as in C11II. All quantities 
expressed in p.u. on a 500 KV, 5000 MVA base. 


0.345 
314 

Z, Excitation System 

= 50 = 0.01 sec. 

3. SVC Voltage Controller 

Kg = 20 Tg = 0.02 sec. 

4. Network Parameters 


1 . Generator 


T* = 6.4 

do 


X = 1.64 

q 

H = 4.2 


X. = 


Ci = 


i) Step up transformer = 0.15 

ii) Transmission (500 KV, 100 km, line constants 
circuit) 


f 1 ■ 
line 

^line 

®line 


0.0135 

0.452 

0.02294 


6. Infinite bus voltagei 

e. Kg ■ O. 


1.0 



p.u. 


are 


per 



APPENDIX C 

POLE ASSIGNMENT WITH OUTPUT FEEDBACK 1 231 


The technique given here is based on the work of Munro and 
Hirbod L363. The method presented in C233 provides a technique 
for the design of full rank compensators for multi variable 
linear systems. Here the technique is restricted to single input 
linear systems. Ue have. 


X = A X + b u 


CC.1) 


and 


y = C X 


(C.E) 


where, x f r”, u f and y f R*^ are respectively vectors of 
state, input and output variables. The r X p open loop transfer 
function matrix G(s) is given by. 


G(s) = C{sl-A) ^ b = C1/AJ <s>3 

- Q 


where 


/X Ti, n-1 

fj (s) = s + at. s + . 

o 1 


N^(s) 


N (s) 
r 


at 


<C.3) 


<C.4) 


and 


Ni(s) = s" + /J.gs" ^ + 


in 


<C.5) 


p^(s) is the characteristic polynomial of GCsX. Consider output 


feedback law 



u(s) 


H(s>. Y(s) 


= u (s) - 

r 


where represents the reference input. 

The problem can be defined as the determination of p X r 
dynamic feedbacka compensator H(s), such that the closed loop 
system, 

H(s) = C1 + G(s) H(s) G(s) (C.7) 

has a desired set of eigenvalues. The order compensator 

H(s) has the form 


H(s) = LA/yi (s)TCM.(s), M„(s) M <s)D 

o 12 r 

where 

/ % k ^ k-l ^ k-2^ 

AJ^(s) = s + s ■*■^2 ® * ^Va 

and 


k k“" 

M. (s) = e.-s + d. .5 
1 lO i1 


+ Q 


ik 


(C.8) 


(C.9) 


(C.10) 


^ (s) is the characteristic polynomial of compensator transfer 
function H( s ) . 

The resulting closed loop system characteristic polynomial a^^Cs), 
defined as: 




n+k ^ . n+k-1 ^ 

= s + d. s + 

1 


. + d 


n+k 


CC.11) 


can be written as [I23D 

r 

AJ .(s) = A^ (s) . p <s) + y N-Cs). M.(s) (C.1E) 

d o c /a 1 I 

i=1 

Complete Pole Placement 

The problem of pole assignment can be defined as given 
p^Cs), N^(s), ... N^(s) and a desired set of closed loop poles, 
find (s), M^Cs) , . . . .M^Cs) of lowest degree which satisfy eqn. 



(C.12). Equating coefficients of like powers of s in eqn. 
(C.1E), we getr 

CX|^3 . ^ (C.13> 

where 
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= C(d.-a 


1>l 




>l 


I (d -a )| Id 
• n n ' ' 


n+1 ' 


n+k 




Equation (€.13) is a set of (n+k) equations in C(k+p)r+k3 
unknown parameters of H(s). The difference vector contains 
the coefficients of the polynomial 


/J .(s) = P^(s) - P„(s) 

d s o 


k 

s 


(C.14) 


A necessary and sufficient condition for the existence of 
solution of eqn. (C.13) is 



Rank CX. D 
k 


(C.15) 


= Rank C 6^ J 
and a solution of eqn. (C.13) is 




(C.16) 


where X®^ is the generalized inverse of Xj^ such that 
CX^D CXI^I CX^3 = CX^D 


Partial Pole Placement 

Consider the case when only q poles are to be assigned, 
t = (n+k— q) poles are allowed to assume arbitrary locations. For 
this case the closed loop systema characteristic polynomial 
has the form 


Pd<s) = 

(s’ + d,,"-' + . 
1 

. . + 

d )<s^ + 

q 

t- 

^ .. 

. .+e^) 








(C.17) 


where d^ , dg 

,....d^ are specified 

and e^. 

Gg r ■ • 

•■n 

are to 

be 

determined. 

Here the difference 

vector <5, 

— k 

can 

be 

obtained 

by 

equating the 

coefficients of 

like 

powers of 

s in 

eqn a 

(C.17) 

and 

{C.12) as 

« 

e . 
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a • a • 



(C.18) 
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where the 

« 

vectors and 


(i = 1 

F ■ ■ 

a a jr 

t) contain 

respectively 

the coefficients 

of 

the polynomials 





Pd ts) = 


. s’- - p^(s) 


k 

s 


p (s) 

q 


(C. i9) 



f 


f 


p"(s) = p (s) . (i = 1, ,t) (C.20) 

Eqn. (C.13) in this case will take the form 


where 

CD3 = C«5^, 6 ^, 6^ □ 

and 

_ r 

e “* i.e^ f e^ V ■ ■ * • » e^ j 
Eqn. (C.21) can be rearranged as 
P* = 

where 

Cxj^D = CX^, -DD 

and 

t t ^ 

Pk = 3 


(C.21) 


(C.22) 


(C.23) 


Eqn. (C.22) is a set of (n-fk) simultaneous algebraic equations 

9 

in Cr(k+p) + k+t3 unknowns. Vector Pj^ contains parameters of 
H(s) and the coefficients of polynomial of unassigned closed loop 
poles. The necessary and sufficient condition for the existence 
of solution of eqn. (C.22) is given in eqn. (C.15) with Xj^, Pj^ 

V V « 

and replaced by Xj^, Pj^ and respectively. 


